Anomalous condensate fluctuations in strongly interacting superfluids 
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We show that the condensate occupation of a superfluid 
Bose liquid quite generally exhibits anomalously large fluctu- 
ations at finite temperatures T 7^ 0. In three dimensions, the 
variance (SNq) of the number iVo of particles in the conden- 
sate scales nonlinearly with volume V like r^y*/'^ at low T, 
generalizing the result obtained by Giorgini, Pitaevskii, and 
Stringari for a weakly interacting Bose gas. In two dimensions 
there is only a quasicondensate whose fluctuations are of the 
same order as the mean value. 
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The experimental realization of Bose Einstein conden- 
sation (BEC) in alkali vapours has stimulated extensive 
theoretical research on the physics of BEC ||l| . It is well 
known that the standard grand canonical ensemble used 
in textbooks gives unphysically large condensate fluctu- 
ations, {5N§) ~ (iVo)2 oc V^, even at T = 0. Much 
recent theoretical work has thus been devoted to the 
investigation of fluctuations in different statistical en- 
sembles |^"0l- With a fixed total particle number, the 
pathologies of the grand canonical ensemble are removed, 
however the condensate fluctuations remain anomalously 



in harmonic traps, where the condensate fluctuations are 
crucial for determining the linewidth of an atom laser 
|0,E2| as recently shown by Graham |1^ . 

We restrict ourselves to Bose liquids in a fixed volume 
Vd — L'^ with periodic boundary conditions. As shown 
by Feynman, the spectrum of excitations as fc — > is 
exhausted by phononlike modes with ojk = cfc. For low 
enough temperatures such that mc^ 3> ksT , only phase 
fluctuations are relevant |Q and the field operator may 
be written as ffsl 
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with no being the bare condensate density at T = which 
would be observed in a small sample, i.e. in the absence 
of low energy excitations. The operator of phase fluctu- 
ations is 
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where m is the particle mass, c the actual velocity of 

sound, and n the mean particle density. The operators 

,t 



1 /rA"r9\ m2TM/^ i 1 i r -inn ■ cl and c.- arc phouou creation and annihilation opera- 

large, {5Nq) (X T^V^''^, at least for an ideal Bose gas m k • V u • • j- • 



a box. Exploiting the fact that BEC in a noninteract- 
ing gas is described exactly by the spherical model, this 
was first realized by Fierz 0. Surprisingly, the same 
behaviour was found to persist in the presence of weak 
interactions as recently shown by Giorgini, Pitaevski and 
Stringari using the Bogoliuobov theory H . 

In this paper, we demonstrate that the anomalous con- 
densate fluctuations are a general feature of any super- 
fluid Bose system, valid for arbitrary interactions. The 
large fluctuations at finite temperature will be shown 
to be a direct consequence of Bogoliubov's well known 
1/fc^-theorem for the static susceptibility Xv>ipik) of su- 
perfluids H]. They are thus related to the existence of 
gapless modes and correspondingly are absent in artifi- 
cial models with a finite gap in the spectrum as intro- 
duced recently in this context 0] . We give a very simple 
derivation similar to the one employed to prove that the 
thermal depletion of the condensate in dimension d ~ 3 
is oc T^ [|0|. For a Bose liquid in d = 2 there is only 
a quasicondensate whose root mean square fluctuations 
{SNq)^/"^ are of the same order than the mean value (Nq), 
both scaling with volume with a temperature dependent 
exponent 1 — r](T)/2. Being a rather general property 
of a Bose condensed system, the anomalous scaling with 
volume or particle number will also be present for BEC's 



tors, respectively. The prime indicates summation over 
all A: ^ with |A:| < A, A being a momentum cutoff (the 
term fc = is omitted because it corresponds to an irrel- 
evant global phase factor which can always be absorbed 
in the bare condensate density no and drops out in any 
gauge invariant quantity as calculated below). The ac- 
tual, renormalized condensate density is obtained from 
the off diagonal elements of the one particle density ma- 
trix, which gives 
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as long as phonon interactions and higher en- 
ergy excitations are neglected. For large systems, 
lim|j|^oo(i^(x)(^(0)) is always small compared to {(p^(x)), 
which is independent of x. Separating off the quantum 
fluctuations of (^ at T = 0, the mean condensate density 
at temperature T may be written as 

no(T) =no(r = 0)e-(<*'(°»lT-<^'(°)H.^o). (4) 

In a strongly interacting Bose liquid, the condensate den- 
sity at zero temperature may be much smaller than the 
total density n, as is the case in superfluid ^He, where 
?io(r = 0)/n~0.1 0. 



In d = 3, the exponent in (^ is independent of the 
system size and proportional to T^, thus describing the 
well known thermal depletion of the condensate flG] 
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In d = 2, the thermal contribution to phase fluctua- 
tions diverges logarithmically with system size |l7l 
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Hence, the condensate density at finite temperature de- 
pends on the system size via no{T) — no(0) {AL)^'^^'^\ 
with ri(T) oc T. In particular, for any finite T the con- 
densate density vanishes in the thermodynamic limit, re- 
flecting the fact that there is no true condensate in d = 2 
in agreement with the rigorous proof given by Hohen- 
berg |l|. 

The operator of fluctuations around the mean field is 
given by 
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In terms of 5'^{x), the number of particles out of the 
condensate is 
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Since the total number of particles is fixed, the fiuc- 
tuations of the condensate occupation Nq may be cal- 
culated from (SNq) = (SN^^f.), i.e. the condensate and 
noncondensate particles act as effectively infinite particle 
reservoirs for each other |g,^ . Using ^ one obtains 



m^) 



2nge-2(^^(0)) 



dxdx' {(p{x)ip{x'))^ 



(9) 



since {(f{x)(p{x')) is small at large distances both in d = 
2 and in d = 3, provided that the temperature is low. 
Evaluating the correlation function at T = 0, we find 
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The fiuctuations are thus normal, scaling essentially lin- 
early with the volume of the system. 

At finite temperatures, it is convenient to express the 
integral in (0) in terms of the static susceptibility Xvvi^)- 
Using the classical version of the fluctuation dissipation 



theorem, which is always applicable at low k because 
fi^k '^ ksT, we find 
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Since the superfluid density n^ which usually appears in 
X^p^p iql coincides with the full density n at zero temper- 
ature for any translation invariant superfluid, we have 
Xff^}^) ~ iT^/nh'^k'^ for A: ^ 0, independent of the inter- 
action strength. Hence, at low temperatures 
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For the weakly interacting Bose gas there is no depletion 
of the condensate at T = to lowest order and thus 
with no{0)/n replaced by one, our result reduces to the 
one given in 0. The somewhat surprising fact that the 
velocity of sound does not appear in (12), may be traced 
back to the well known low fc-divergence p5lu9| 
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of the momentum distribution at low but flnite tempera- 
tures which, apart from the ratio nQ{Q)/n, only depends 
on the bare single particle energy e^ — {hk)'^/2m. In 
order to make this connection more explicit, we note 
that the large distance behaviour of the phase correla- 
tion function which appears in (9) is directly related to 
the momentum distribution nk of the noncondensed par- 
ticles at low k via [ |l5[ 
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This shows that equ. (9) can be written in the simple 
form 
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which immediately gives (12) using the result (13) for 
the momentum distribution at finite temperature. The 
anomalous scaling of the condensate fluctuations with 
volume is therefore a direct consequence of the 1/fc^- 
divergence of the momentum distribution or - equiva- 
lently - the related slow l/r-decay of the one particle 
density matrix at large distances to its limiting nonzero 
value no(T). Incidentally this argument also explains 
why the result for an ideal Bose gas, where rio(O) = n 
trivially, is just twice the result (12). Indeed the mo- 
mentum distribution for noninteracting Bosons at small 
nonzero k is nj, —> kBT /sk (the factor two difference 
compared with (13) is due to the fact that at low k the 
quasiparticles of the interacting system are just an equal 



weight superposition of a bare particle with momentum 
k plus a bare hole with momentum — fc |20[|). Using the 
standard result for the occupation number fluctuations 
of ideal Bosons we have 
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To leading order the fluctuations in the interacting and 
the nonintcracting case are thus simply related by a fac- 
tor 2(7io(0)/2n) . Contrary to the argument in m, the 
similar behaviour of the fluctuations in an ideal and in- 
teracting Bose gas is therefore not accidental. 

As implied by (|l^), at T 7^ the fluctuations of the 
ground state occupation are anomalous for any interac- 
tion. Expressing the result in terms of the standard ther- 
mal wavelength At = h{2T:mkBT)^^^^ we have 
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with a universal constant B = ^- '|r?|^'*/27r^ = 0.8375 
whose numerical value is readily obtained from the lat- 
tice sum ( [l2|) (this result is roughly a factor five smaller 
than the corresponding one given in ref. M). Note that 
contributions k ~ 2T:n/L with one or two of the compo- 
nents of n G Z^ being zero should be included in the sum 
in the case of periodic boundary conditions as assumed 
here. For Dirichlet boundary conditions in a box with 
hard walls the corresponding constant in an ideal Bose 
gas is in fact smaller by a factor 0.864 pi] ]. 

In d = 2, the one particle density matrix decays to 
zero like l/r^^'^'> at any finite temperature and the asso- 
ciated momentum distribution Uk is thus proportional to 
^-(2-17) at small k. The resulting scaling of the conden- 
sate fluctuations with system size is temperature depen- 
dent 
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increasing like the square of the average condensate num- 
ber (SNq) oc (No)^- The fluctuations in d = 2 are thus 
much larger than those in d = 3. The fact that the rms 
fluctuations of the condensate are of the same order as 
the average is yet another indication that the condensate 
in d = 2 is not well defined. 

In conclusion, we have shown that the anomalous fluc- 
tuations of the condensate occupation at finite temper- 
ature, found previously for both ideal and weakly in- 
teracting Bose gases, are a general property of any su- 
perfluid with arbitrary interactions. . The anomalous 
behaviour of these fluctuations may be traced back via 
Bogoliubov's l//c^-theorem to the corresponding low k- 
singularity in the momentum distribution. As has been 
shown in H, the condensate fluctuations in harmonic 



traps are similarly determined by the associated low lying 
collective excitations. As a result, the anomalous scaling 



iSN^] 



rp2j\^4/3 -^j^i^ ^i^g particle number N is also 



expected to be valid quite generally 
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